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Abstract 

Accurate non-perturbative calculations of glueballs are performed using light-front quantised 
SU{N) gauge theory, to leading order of the expansion. Based on early work of Bardeen 
and Pearson, disordered gauge-covariant link variables M on a coarse transverse lattice are used 
to approximate the physical gauge degrees of freedom. Simple energetics imply that, at lattice 
spacings of order the inverse QCD scale, the effective light-front Hamiltonian can be expanded 
in gauge- invariant powers of M: a colour-dielectric expansion. This leads to a self-consistent 
constituent structure of boundstates. We fix the couplings of this expansion by optimising 
Lorentz covariance of low-energy eigenfunctions. To lowest non-trivial order of the expansion, 
we have found a one-parameter trajectory of couplings that enhances Lorentz covariance. On 
this trajectory the masses of nearly-covariant glueball states exhibit approximate scaling, having 
values consistent with large- A'' extrapolations of continuum results from other methods. There 
is very little variation with N in pure Yang-Mills theory: the lightest glueball mass changes by 
only a few percent between SU{3) and 5*17(00). The corresponding light-front wavefunctions 
show an unconventional structure. We also examine restoration of rotational invariance in the 
heavy-source potential. 
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1 Introduction 



There are few, if any, efficient methods for tackUng relativistic strongly-bound states in generic four- 
dimensional gauge theories. The canonical example is QCD, where non-perturbative theoretical 
calculations have rarely been ahead of experiment. Future progress in particle physics is likely to 
hinge upon a detailed theoretical understanding of these questions. This has led some theorists to 
develop Hamiltonian quantisation on a hght-front (l], g. In the presence of suitable high-energy 
cut-offs, the light-front vacuum state is trivial and wavefunctions built upon it are Lorentz boost- 
invariant. 

In particular, Brodsky and Lepage |^ and Pauli and Brodsky Q have urged the development 
of the light-front quantisation of QCD (LFQCD). More recently, Wilson et al. ||] have clarified the 
physical principles underlying LFQCD and suggested a weak-coupling calculational framework. In 
this paper we develop an alternative framework which appears promising: a light-front quantisation 
of lattice gauge theory . Calculations that we have already performed with this method ^, |o| , 
for non-Abelian gauge theories in 2 -I- 1 dimensions, were surprisingly successful in comparison to 
results from traditional Euclidean lattice path integral simulations (ELMC) They have also 

produced new results in the form of the light-front wavefunctions, the starting point for investigation 
of virtually any physically interesting observable. Encouraged by this success, we investigate here 
the glueballs and heavy-source potential in 3 + 1-dimensional gauge theory without fermions. 

The well-known triviality of the light-front vacuum may be reconciled with the conventional 
picture of a complicated QCD vacuum. In light-front co-ordinates, vacuum structure is carried 
by an isolated set of (infinitely) high energy modes, which are removed by the cut-off. According 
to standard lore, one would expect that appropriate renormalisation of the Hamiltonian and other 
observables would recover the information excised by the high-energy cut-off. In this way, effects nor- 
mally associated with the vacuum, such as spontaneous symmetry breaking, must appear explicitly 
in the Hamiltonian. But the non-perturbative renormalisation group (RG) formalism to systemati- 
cally compute the necessary counter-terms does not yet exist. Efforts are being made to formulate 
a perturbative light-front RG for QCD JT^ , but the full vacuum structure must necessarily appear 
via counter-terms that are non-analytic in the perturbative couplings. Faced with the problem of 
finding non-perturbatively renormalised Hamiltonians, we will use symmetry as our guide. 

The first step is to choose the most general set of Hamiltonians which respect symmetries of the 
theory that are not violated by the cut-offs. After truncating this set according to some reasonable 
criteria, we non-perturbatively test the low-energy eigenfunctions and eigenvalues for restoration of 
the symmetries violated by the cut-offs. These tests, rather than RG transformations, are used to 
explore the space of Hamiltonians. Since we maintain gauge invariance, the symmetries in question 
involve only Lorentz covariance. We will remove all but one cut-off and find that our truncated 
space of Hamiltonians contains a unique one-parameter trajectory Ts on which Lorentz covariance is 
greatly enhanced. The validity of Ts is confirmed by the approximate scaling of low-energy physical 
quantities along this trajectory. Ts will be used as the basis for extracting cut-off independent 
results. This is a first-principles approach, since no data are taken from experiment (aside from the 
overall scale in QCD). 

Apart from Lorentz invariance, we must address the issue of gauge invariance. One knows that 
gauge invariance can be maintained with lattice cut-offs and compact variables [llSl llj]. However, in 
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light-front coordinates a high-energy lattice cut-off can only be applied to the transverse directions; 
the null direction on the initial surface must remain continuous. Also, if one were to employ the usual 
formulation of lattice gauge theory with degrees of freedom in SU{N), it is not straightforward to 
identify the independent degrees of freedom that are essential for canonical Hamiltonian quantisation. 
The tricks of 'equal-time' Hamiltonian lattice gauge theory in temporal gauge do not carry over 
to light-front Hamiltonian lattice gauge theory in any convenient way p^ . One docs not have to 
choose lattice variables in SU{N) however. It was noted long ago by Bardeen and Pearson that 
lattice variables M in the space of all complex N x N matrices were physically more appropriate on 
a coarse lattice. Gauge invariance is maintained and it is straightforward to identify the independent 
degrees of freedom. The penalty is that one is too far from the continuum to use weak-coupling 
perturbation theory. But in the case of light-front Hamiltonians, this may be of limited use anyway. 
Bardeen et al. Q suggested that on coarse lattices one could expand the light-front Hamiltonian in 
powers of M (a kind of strong-coupling expansion). The validity of this expansion is subtle however, 
since it rests on the dynamical properties of light-front Fock states in this kind of theory, in particular 
the weakness of couplings between sectors with different number of partons []l6| . 

Expanding the light-front Hamiltonian of SU (N) gauge theory to leading non-trivial order in M 
and 1/N, we find the trajectory of couplings 7^ mentioned above. In the region of coupling space 
we are able to investigate, the transverse lattice spacing on 7^ is found to be always greater than 
about 0.65 fm. However, Lorentz covariance and scaling are present to sufficient accuracy that we 
can make direct estimates of the continuum values of low-energy observables. Q 

The organisation of the paper is as follows. In the first part of the present work, we give a 
systematic treatment of tranvserse lattice gauge theory. In Section ^ we construct generic light-front 
Hamiltonians on the transverse lattice and introduce approximation schemes for studying them; in 
particular, we review the 'colour-dielectric' expansion in powers of M that leads to a constituent 
picture of boundstatcs. By working to leading order in 1/iV, the transverse dynamics dimcnsionally 
reduce (for a coarse lattice) the problem becoming mathematically equivalent to a 1 + 1- 

dimensional gauge theory with adjoint matter As well as the lattice cut-off, which can only 

be used in transverse spatial directions, we also employ null-plane boundary conditions, specifically 
Discrete Light-Cone Quantisation (DLCQ) and Tamm-Dancoff cut-offs for transverse degrees of 
freedom. The latter gauge-invariant cut-offs can be removed by extrapolation to give finite answers 
for fixed transverse lattice spacing, the only remaining cut-off. 

The renormalisation of the theory at fixed lattice spacing is then accomplished by optimising 
Lorentz covariance, discussed in Section ^. This amounts to a choice of 'metric' in the space of Hamil- 
tonians, based on explicit calculations of observables. Results from our calculations are presented 
in Section ^; we will not describe in any detail the calculational techniques, but refer the reader 
to our previous papers. Some preliminary results from the giueball calculations were presented in 
Ref. 0, but the new material here is more accurate and extensive, including a thorough analysis 
of the space of couplings, the giueball masses and wavefunctions, and the heavy-source system. 

We end this introduction by quoting the most readily understandable of our main results: the 
lowest giueball masses. We estimate the lightest giueball in SU{(X)) gauge theory, the J'^^ — 
to be at a mass Mq++ — 3.50 ± 0.24^/0-, where cr is the string tension. This is essentially 

^As with other 'strong-coupling' approaches, the question of the formal connection to the continuum limit is not 
directly answered. 
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Figure 1: The variation of glueball masses with N (pure glue). ELMC predictions are continuum 
ones for iV = 2, 3 |jl^, |l^, ^ and fixed lattice spacing estimates for iV = 4 The dotted lines 
are to guide the eye and correspond to leading linear dependence on 1/7V^. 



indistinguishable from the result for 5[/(3) pure gauge theory established rigorously by ELMC in 
recent years For the 2++ tensor state, not all components of the multiplet are yet behaving 

covariantly; but from those components which are reliable, we estimate M2++ — 4.97± 0.43a/ct. The 
vector l+~ has Mi+~ = 5.57±0.4-\/ct. We find no fight pseudo-scalar '", but suspect our candidate 
for this has a large error. These results are summarised in Figure |l|. Other recent estimates of glueball 
masses have been obtained in SU (3) Hamiltonian gauge theory |^l| and based on extensions of the 
Maldacena conjecture for large- iV gauge theory though we consider those results to be less 
rigorous than ours. 



2 Transverse lattice Hamiltonians. 

2.1 Energy-momentum 

In 3-)- 1 spacetime dimensions we introduce a square lattice of spacing a in the 'transverse' directions 
X = {x^ , x^} and a continuum in the x^} directions. In light- front (LF) co-ordinates x^ = {x^^ ± 
x^)/y/2, we treat x'^ as canonical time and place anti-periodic boundary conditions on x^ ^ x^ +£. 
Both 1 /a and £ are high-energy cut-offs for the LF Hamiltonian P~ = {P^ — P^) that evolves the 
system in LF time The Lorentz indices n,^ ^ {0, 1, 2, 3} are split into LF indices a, /? G — } 
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and transverse indices r, s G {1,2}. 

The gauge field degrees of freedom below the cut-offs are represented by Hermitian gauge poten- 
tials Aa{x) and complex link variables Mr{x). We also introduce heavy scalar sources 0(x). On the 
transverse lattice, Ao,(x) and (^(x) are associated with a site x, while Mr(x) is associated with the 
link from x to x + af. Each 'site' x is in fact a two-dimensional plane spanned by {a;+, x~}. These 
variables transform under transverse lattice gauge transformations U G SU {N) as 

A„(x) ^ C/(x)A„(x)f/t(x)+i(9„C/(x))C/t(x) 

M^(x) ^ ?7(x)M^(x)?7t(x + ar) (1) 
(A(x) ^ [/(x)</-(x) . (2) 

Since it will be possible to eliminate Aa by partial gauge-fixing, M and represent the physical 
transverse polarisations. The simplest gauge covariant combinations are M, (p, F"^, detM, D"M, 
D"<p, where the covariant derivatives are 

DaMrix) = {da + i^a(x)) Mr(x) - zMr(x)Aa(x + af) (3) 
Dacp = dacp + iAa(t) . (4) 

From these we wish to construct general LF Hamiltonians invariant under gauge transformations 
and those Lorentz transformations unviolated by the cut-offs. 

To proceed, we must make some assumptions about which finite sets of operators to include in 
the calculation. Since symmetries will be tested explicitly, a poor choice of operators would show 
up later on. The following criteria were used to select operators in P~ for pure gauge theory: 

(A) Canonical quadratic form for P+; 

(B) Nai've parity restoration as £ — > oo; 

(C) Transverse locality; 

(D) Expansion in gauge-invariant powers of M. 

Each of these criteria deserves some explanation. The last three can all be straightforwardly checked 
in principle by systematically relaxing the condition. 

(A) Of the generators of the Poincare group {P+, P", P'^, M^^}, the subset {P^ P+, M+'^, M^^, 
M"* } can usually be made kinematic. That is, they can be made independent of interactions, 
quadratic in the fields. The imposition of a lattice cut-off can spoil this property, especially 
if one wants to maintain gauge invariance. We try to maintain as many operators as possible 
in kinematic form, consistent with gauge invariance. Further details are given below, but 
for now we note that condition (A) can be satisfied in light-front gauge A- = by using a 
Lagrangian containing the only ^'''-dependent gauge-covariant term quadratic in link fields. 
DaMr{'x.){D°' Mr{'x.)y . Although higher order terms in M cannot be ruled out, only with a 
quadratic term is quantisation straightforward and, even then, only in the LF gauge A- = 0. 

(B) We will extrapolate to the £ ^ oo limit in the longitudinal direction, deriving P~ from a 
Lagrangian including only dimension 2 operators with respect to {x"*" , x~ } co-ordinates. It has 
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been noted that functions of dimensionless ratios of longitudinal momenta p+ can appear in 
couplings of boost-invariant LF Hamiltonians Q. However, these functions must be strongly 
constrained by LF parity <-*■ x~ . Parity is a dynamical symmetry on the light-front (it 
does not preserve the quantisation surface) and is difhcult to check explicitly. We will assume^ 
that, if is derived from a Lagrangian with naive x+ <-^- x" symmetry (^^-independent 
couplings), then parity is correctly restored in transverse lattice gauge theory in the limit 
C oo. Although DLCQ treatments of analogous 1 -t- 1-dimensional gauge theories have been 
successful under similar assumptions |]l6|, in the present case there are many more 

possible operators that are ruled out by conditions (A) and (B) combined. 

(C) In products of gauge-invariant operators on the transverse lattice, each term of the product 
can be arbitrarily separated in the transverse direction. We assume some transverse locality 
by restricting products of gauge invariant operators to share at least one site x. 

(D) After fixing to LF gauge = and eliminating the resultant constrained field Fock space 
states will consist of link-partons derived from the Fourier expansion of M . For sufficiently 
large a, M remains a massive degree of freedom and there is an energy barrier for the addition 
of a link-parton to a Fock state. Operators of order in will connect Fock states differing 
by at most p — 2 link-partons. By expanding P~ to a given order of M in this regime (the 
colour-dielectric expansion) we therefore cut off interactions between lower-energy few-parton 
states and higher-energy many-parton states. Note that all energy scales are still allowed 
however, since we take the C oo limit, enabling highly virtual small p+ partons to appear. 
The advantage of a cut-off on changes of parton number is that it organises states into a 
constituent hierarchy, consistent with energetics. 

The Lagrangian density satisfying condition (A) can be written 

Lx = :D,M,(x)(D"M,(x))t - 14 - C/x (5) 

where the 'potential' has a purely transverse part and a mixed part ?7x. Up to 4th order in M, 
the purely transverse part is 



Tr|Mi(x)M2(x + ai)M|(x + a2)Af2^(x)| + c.c. 



^ Tr{M,Aft} + ^ ^ (det[Af,] + det[Mt] 



y 4^7 Tr{M,MtM,Mt} + ^ Tr{ A'/,(x)M,(x + af)M}{ii + af)Aft(x)} 



A3 r , , . ,+ -,x2 A._ 

N 

a=±2,a' = ±l 

4A, 



+ E^(Tr{M.Aft})%-^ J2 Tr{A^>/.A//tA/.} 
^2^2Tr{Ma4}Tr{Af2A4} . (6) 



^This assumption is not warrented in the presence of finite-mass fermions. The subtlety lies in p+ = vacuum 
modes that are not recoveredeven as £ — ► 00. For Yukawa interactions of finite-mass fermions, parity certainly isn't 
recovered in the naive way [bsl . 
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The det[M] term can be dropped if > 4 and can always be neglected in the large- iV limit. The 
mixed part, to the approximation we will need, can be written 

C/x - ^ Tr {E^P^^F^^{y)Pl^} - G^f{E^) , (7) 
y 

where E^^ is a non-dynamical pseudoscalar adjoint field at site x, while Pxy is a linear combination 
of Wilson lines in M, each from from x to y (for gauge invariance). Upon integrating out E, the 
function f{E) = Tv{E'^] + 0{E^) gives a simple ^Tv{F°'f^ F^p] term if only the E"^ part is 
retained and we set Pxy = <5xy0 This will be the only term needed for pure gauge theory to 0{M'^) 
once Aa has been eliminated. We gave the more general form in Eqn. (^ since it will be relevant in 
the heavy-source analysis. 

From the above Lagrangian, four of the usual seven kinematic Poincare generators derived canon- 
ically from the energy-momentum tensor T^'' remain gauge invariant and can be made kinematic 
by LF gauge choice yl_ = 0. They are, at x+ = say, 

P+ 2 /"dx- ^ Tr{a_A/^(x)a_Af^(x)t} , (8) 

= 2 / dx" ^ Tr { a_ Ms (x) d- (x) ^ } (9) 
= -2y"da;-^ (x'' ^J""^) Tr{9_A/,(x)9_A/,(x)t} (10) 



M 



x,s 



We will use these to define states of definite momentum {P+,P}. The other three would-be kine- 
matic Poincare generators in light-front formalism, {P^ , Af^^}, are not gauge-invariant when derived 
canonically from T'^'^ because of the lattice cut-off; the simplest gauge-invariant extensions are no 
longer quadratic in fields, even in LF gauge. 

Of the dynamic generators, the most important, and the only one we explicitly treat, is the 
light-front Hamiltonian itself 



} , (11) 



= 1 dx- Y (^x - Tr{A+(x) J+(x)} - Tv{d^A+d^A, 
J+(x) = i Y (A^r(x) d- Aft(x) + M}{x - af) d- Mr{-K ~ af)) . (12) 

r 

Rather than trying to directly construct an approximate realisation of the Poincare algebra at finite 
cut-off, we will minimise cut-off artifacts by optimising restoration of Lorentz covariance in low- 
energy eigenstates of the Hamiltonian p- . 

In the light-front gauge, A+ is a non-dynamical variable and eliminating it introduces non-local 
interactions thus 



P- = 



where /d- = d_^{J^). There is still a residual x -independent gauge invariance generated by 
the charge / dx- . As originally shown in Refs. ^, |7j, finite energy states \^) are subject to 

•^The case -fxy — ^xy is related to the generalised 2DQCD of M. Douglas e£ al. |26| 
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the gauge singlet condition J dx~ J~^\'i') = 0. In the large- limit, this means that Fock space at 
fixed is formed by connected closed loops of link variables M on the transverse lattice (the x~ 
co-ordinate of each M is unrestricted). 

2.2 Quantisation 

The dynamical problem is now to diagonalise at fixed total momenta {P^,P}. A convenient 
basis consists of free link-partons obtained from the Fourier modes a(fc+.x) of AI in the x^ co- 
ordinate 

M,(x+ = 0,a;-,x) = -= / ^= ( a_,(fc+, x) e"''^^" + 4(fc+, x) e*'=^ ) . (14) 



VK^ ^ ' 

The Fock space formed from creation operators in the large- limit, and other details of the 
calculation in this Fock space, including construction of states of definite momentum, have been 
described elsewhere [|[ |l^. We applied both DLCQ and Tamm-Dancoff cut-offs in Fock space, 
extrapolating both of these at fixed values of the couplings in Eqn. (||).0 

Low-energy eigenfunctions of P~ {id est glueballs) are to be tested for Lorentz covariance, the 
couplings appearing in Eqn. (^) being tuned to minimise covariance violations. Since G'^N, with 
dimension (energy)^, is consistent with 't Hooft's large- iV limit |2^, we will use it to set the dimen- 
sionful scale. Thus, the dispersion relation of a glueball can be written 

2P+P- = G^N (Ml + Mla^\Pf + 2Mla^P^P^ + 0{a^\Pf)) . (15) 



For each glueball, Ai^ are dimensionless functions of the couplings which, for a given set of couplings, 
are extracted by expanding eigenvalues of P^ in aP. A truly relativistic state must have an isotropic 
speed of light 

a^G^NMl = cl^ = 1 
a^G^N{Ml+Ml) = cls = 1. (16) 

Con is the speed of light in direction x = (1,0); Cos is the speed of light in the direction x = (1,1). In 
addition, the an-harmonic corrections at 0(a'*|P|^) must vanish, but we will not use this condition 
directly. 

The ratio of transverse to longitudinal scales is set by the dimensionless combination a^G^N . 
This can be deduced from two measurements of the string tension ct, for example |2|]. The first 
measures the mass squared of winding modes on a periodic transverse lattice, fitting the groundstate 
to n^a?(j^ for winding number n. The second measures the heavy-source potential in the continuous 
x^ direction, fitting the groundstate to ctj^R for separation i?.^ Demanding ctt^ = (t^ = cr and 
expressing eigenvalues in units of G'^N, we may eliminate a and deduce a^G^N. We may also 
deduce G^N , and hence all dimensionful quantities, in units of a. Some refinements we have made 
to the calculation of the heavy-source potential, since our previous work, are described in the next 
subsection. 



*The largest cut-offs used were a K = 26/2 DLCQ resolution and 8-link Tamm-DancofT cut-off in the glueball 
sector, while a K = 70/2 and 4-link cut-off was used in the heavy-source sector. The latter sector has an additional 
cut-off Pmaxi discussed in Ref. |lo| , which was varied up to 7. 

^One could also use the heavy-source potential in transverse directions to find aT, but it is more accurate to use 
winding modes for this. 
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2.3 Heavy-source potential Vqq 

The non-relativistic heavy-source potential Vqq provides us with a number of important pieces of 
information. The string tension a wiU be used as our basic physical unit for dimensionful quantities. 
The deviations of the potential from rotational invariance will help to fix couplings in the theory. 
We find that, in general, both relativistic invariance in the glueball sector and rotational invariance 
in the heavy-source sector are needed to accurately pin down the effective LF Hamiltonian. 

The heavy-source sector was first studied on the transverse lattice in 2-1- 1 dimensions by Burkardt 
and Klindworth [ p8| . They used the most rudimentary approximation to the Lagrangian, quadratic 
in the link and heavy-source fields. Along both the continuous and lattice spatial directions linear 
potentials arise, which can be given the same slope a by adjusting the link-field mass /z^ for given 
a/y/a. The authors of Ref. observed that the contours of the potential in the spatial plane 
were 'square' a.t fi^ ^ oo {a ^ oo) and 'circular' at /i^ = (a ^ 1/y/a). This rounding of the 
potential indicates that rotational invariance is improved as the lattice spacing is reduced, as one 
would hope. However, there is still a potentially large violation of rotational invariance due to the 
fact that any discrepancy between the tensions measured in continuum and lattice directions can be 
masked by adjusting a. Only by examining another sector where the lattice spacing value is used, 
such as covariance of glueballs, can one tell whether the potential is truly circular or merely oval. 
Using more accurate higher order LF Hamiltonians in the pure glue sector [|lO[ we found that 
the values of the lattice spacing deduced from rotational invariance of the heavy-source potential 
were not very compatible with those needed for relativistic covariance of glueballs. Therefore one 
needs to also improve the accuracy in the heavy-source sector. We take some steps in this direction 
in the following. 

Our light-front treatment of heavy-source fields (f>i , which we take to be scalars in the fundamental 
representation, has been introduced in Ref. We will work with a LF Hamiltonian in LF gauge 
A- = to order in the pure glue operators, order in the heavy-source operators, leading 
order of the heavy-source mass p, together with the conditions discussed in Section ^. This LF 
Hamiltonian for heavy sources with non-zero spatial separation can be derived from the Lagrangian 
density 

ix - ^ Tr|:D„Af,(x) (:D"Af,(x)) ^1 - ^ Tt{F"^F„p} -V^ + {D^<f>)^ D^<jy - p^^t^ 

Trl^^"'^ F^p Wr] - ^ Tr{Mt(x) F"/5(x) M.(x) F^p{^ + ar)] , (17) 

where 

Wr = {MlMr + MrMl) . (18) 

In addition, there are combinations which are suppressed in the p ^ oo limit, such as (f)^Wr4>, 
(j)^ F"^ Fai3(j), which we did not include in this work.^ 

After gauge fixing A- — 0, eliminating A+ in powers of M , and discarding the higher orders in 
M, the LF Hamiltonian resulting from ( p^ is 




^We investigated the first of these terms in ear Her work by scaling its coefficient with p so that it survived the 
p — > oo fimit. It had Uttle afi'ect on the string tension but seemed to help short-distance rotational invariance. 
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with 



J+t ^J++i<t> d- 0t (20) 

The diagonalisation of this Hamiltonian in Fock space proceeds in a similar way to the gluebaU 
sector (see Ref. 
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The selection (|l7l ) is more general than was used in previous calculations. Physically speaking, 
the new terms ti and T2 generate oscillations of the flux string perpendicular to x^-separated sources 
(via enhanced pair production of links) as required of a rough gauge string . This partially screens 
the linear potential in the x^-direction to produce a rotationally invariant a. However, while we find 
that these new terms improve the interaction between heavy sources at zero and one-link transverse 
separation, small- a;^ separations and wider transverse separations are left unimproved. This will 
result in a loss of off-axis rotational invariance, which could only be remedied by including terms 
at O(M^) et cetera. The underlying difficulty in getting rotational invariance may be related to 
the existence of a roughening transition separating the large a regime from the continuum in the 
heavy-source sector. 



3 Space of Hamiltonians 

Our procedure for extracting Lorentz-covariant cut-off independent results is non-standard, so we 
take some time in this section to describe the procedure we used for the present calculation. 

3.1 Topology 

Abstractly, we have a set of cut-off Hamiltonians containing in principle an infinite number of 
operators Oi with couplings gi. For simplicity, we consider only one high-energy cut-off 1/a, all 
others having been removed by extrapolation.^ Ti. is the infinite-dimensional space of the gi, the 
space of Hamiltonians. 

One could in Ti define renormalisation group trajectories as trajectories along which, say, the 
eigenvalues of the Hamiltonian are invariant. We are particularly interested in a trajectory on 
which Lorentz covariance is restored — the 'Lorentz trajectory' T — for observables containing 
no explicit momenta above the cut-off. Perry and Wilson |^ have emphasised the idea, albeit 
in weak-coupling perturbation theory, that when a symmetry is restored in a cut-off theory, the 
number of couplings running independently with the cut-off is reduced. For Lorentz symmetry in 
pure gauge theory, for example, this coherent running of the couplings should pick out T C 7i. 
A basic hypothesis of our approach to pure SU (N) gauge theory on the transverse lattice will be 
that a single one-parameter Lorentz trajectory extends from the continuum to quite low cut-offs. In 
order to test the hypothesis, we will specify conditions for an approximate Lorentz trajectory 7^ in a 
finite-dimensional subspace Hs C H. The subspace Hs arises from the conditions we placed on the 
Hamiltonian in Section ^. The most straightforward way to define Tg is in terms of a finite selection 
of observables contributing to a t^st for deviations from Lorentz covariance. At the same time, 

^This extrapolation in itself may require tuning of the gi in general. However, such tuning is not necessary for 
pure gauge theories within the approximations we make. 
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Figure 2: The space Ti. = {51, 52, 53, ■ ■ ■} of Hamiltonians and its subspace Hs = {51,52}- T is 
the exact Lorentz trajectory, 7^ the approximation to it in Tig ■ C is the intersection of T with the 
continuum 1/a = 00. The dashed hne is a possible phase boundary in Tig separating 71 from the 
continuum. 

we beheve that it should be possible to develop tests of Lorentz covariance that do not depend on a 
specific choice of observables, exempli gratia ensuring closure of the Poincare algebra. But we have 
yet to carry out this development. 

Figure ^ illustrates the space of Hamiltonians we have in mind. The space Ti. = {gi, 32, 53, • • •} 
represents the infinite-dimensional space of couplings. The cut-off 1/a decreases along the Lorentz 
trajectory T as one moves away from the intersection of T with the continuum, C. A finite- 
dimensional subspace Tis is denoted by the 53 = plane. T will not lie entirely in Tig in general, 
so the best one can do is investigate approximations in Hs- We search Hs over a certain range of 
cutoffs for a trajectory 7^ which minimises the distance to T, id est minimises violations of Lorentz 
covariance. It is not sufficient to simply find an isolated point in Tig at which violations of Lorentz 
covariance are minimised. This may have little to do with T. One must establish an entire trajec- 
tory. If T does not exist for the range of cut-off investigated, or lies far from Hs, one will be unable 
to find an unambiguous trajectory 7^, and the method fails. 

The metric to determine the distance from T is the test for Lorentz covariance, which involves 
both a choice of observables and their relative weights. Even if T exists, a random choice of 
observables and their weights will make it difficult to identify. Consistent with the colour-dielectric 
expansion, we gave higher weights to low energy observables dominated by a few link-partons and 
made sure that the number of observables was always ^ dim [Hs] ■ 

Having convinced oneself that 7^ is uniquely established for the metric employed and the range 
of cut-off investigated, one can use it to estimate the value of observables on T. Observables will 
exhibit approximate scaling along 7^, by virtue of its nearness to T (which is an exact scaling 



11 



trajectory). The values of these observables will be close to the values on T, and the violations of 
scaling and covariance as one moves along 7^ can be used to estimate the systematic error from use 
of a finite-dimensional space of Hamiltonians. As one enlarges the space Tig, exempli gratia by going 
to higher orders of the colour-dielectric expansion, one can get closer to T and the errors reduce. 

This procedure is to be contrasted with the traditional perturbative 'improvement' program of 
introducing irrelevant operators in a cut-off theory. There, one derives the couplings via perturbative 
renormalisation group transformations, solves the theory, and extrapolates to the continuum. This 
is the approach to LF Hamiltonians being pursued by a number of authors [^2j. Instead, we use sym- 
metry directly to construct coupling trajectories, making the assumption that symmetry uniquely 
fixes the values of observables in QCD, once the overall scale has been specified. In practice, we do 
not construct renormalisation group trajectories and do not extrapolate to the continuum. 

We do not try to extrapolate the results to a = because in Hs there may be barriers to 
the continuum (exempli gratia roughening or large- TV transitions). Our use of disordered massive 
link variables Af , together with the concomitant expansion of in powers of M, means that the 
transverse lattice spacing a is quite coarse in the region of Ti. we search. However, by definition, 
there can be no phase boundaries separating T <Z Ti from the continuum. Provided a suitable T 
exists, and we can get close enough to it, the interpolation of phase boundaries only breaks the 
connection between scaling violations in each phase of Tig , not the estimate of the observable itself. 

3.2 Topography 

We now construct a set of topographical charts of Hs that help one to decide whether a unique 
Lorentz trajectory exists. It is convenient to use G^N to set the dimensionful scale wherever possible, 
and introduce dimensionless versions of the other couplings as the co-ordinates of Tig 

2 _ /^^ , _ At _ Tj , _ [3 ,„ . 

We employed two ways to search Hg. 

Method 1 For a given m^, this minimises over all couplings {^i, ^2, ^3, ^4, ^5, to give a one- 
parameter trajectory in Tig. This amounts to a direct determination of Tg. 

Method 2 For a given m and gi £ {Ii,h,l3,l4,l5,b}j this minimises over all the remaining 
couplings. There is a chart for each gi which can be displayed as a contour plot with heights 
xinin- The existence of T should show up as a unique valley on each chart, universal with 
respect to small changes in the form of the metric, the bottom of the valley coinciding with 
the trajectory Tg found by Method 1. 

With Method 1, in a reasonable amount of time we could extrapolate the DLCQ and Tamm-Dancoff 
cut-offs before using a routine to iteratively search for the minimum value of x^, setting the mass scale 
from first principles via string tension measurements. Because of its more time-consuming nature, 
a number of compromises were necessary to obtain decent results by Method 2: the DLCQ and 
Tamm-Dancoff cut-offs were fLxed;^ to avoid the extra complication introduced by the heavy-source 
analysis, the mass scale was fixed by setting Ai{0~^^) = 3.5y/a, a result which had been deduced from 

= 6/2 and a 6-link truncation. 
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m 


h 


I 


2 


h 




k 


h 


b 


h 


h 


y 

A. 


2 


0.0444 


0.012 


-0 


068 


437 





097 


182 


0.429 


-0.795 


-0.702 


13 


36 


0.0890 


-0.005 


-0 


075 


426 





118 


197 


0.504 


-0.908 


-0.864 


14 


28 


0.1342 


-0.035 


-0 


079 


463 





132 


211 


0.575 


-0.999 


-1.008 


15 


99 


0.1803 


-0.038 


-0 


090 


374 





158 


143 


0.688 


-1,086 


-1.094 


18 


76 


0.2275 


-0.037 


-0 


101 


62.4 





212 


18.8 


0.858 


-1.205 


-1.189 


21 


14 


0.2765 


-0.048 


-0 


137 


111 





208 


22.4 


0.962 


-1.147 


-1.185 


25 


28 


0.3275 


-0.007 


-0 


144 


60.7 





163 


9.6 


1.010 


-0.788 


-0.734 


31 


87 


0.3812 


-0.007 


-0 


195 


12.3 





308 


2.28 


1.381 


-1.207 


-1.302 


31 


83 


0.4384 


-0.047 


-0 


175 


50 





331 


12.5 


1.530 


0.588 


-0.635 


43 


35 


0.500 


0.319 


-0 


313 


822 





367 


318 


1.906 


-1.42 


-1.605 


59 


78 



Table 1: The trajectory which minimises the test of Lorentz covariance. 



Method 1 on 7^ (see later); we used fixed values ^3 = 100 and Z5 = 10, again motivated by the typical 
values on 7^ from Method 1. The very large couplings ^3 and ^5 have a special significance which 
is explained in Rcf. ]lO[ |. They can be set to virtually any very large value for a Lorentz-covariant 
theory. 

The first step in the process requires some experimentation with different variables. In 
general, we used the isotropy of glueball dispersion relations, rotational invariancc of the heavy- 
source potential, and Lorentz multiplet structure of the low-lying glueball spectrum as the ingredients 
in x^. After some fine-tuning of weights, we found a very clear signal for a trajectory 7^ running 
from small to large m on which the was significantly reduced to about one per effective degree 
of freedom. Fig. ^ and Table |^ show the trajectory found by Method 1. Further details of the 

choice are available in a data file 40 variables were used in x^, though 15 of these were 

given significantly more weight than the others, corresponding to groundstate observables in each 
symmetry sector. The lowest values of x^, in the range 10-20, occur on a smooth trajectory through 
the lowest four non-zero values of m sampled, above which 7^ begins to break up somewhat. Fig. ^ 
shows the global topography of coupling space via Method 2, confirming that the trajectory found 
in Method 1 sits at the bottom of a well-defined, unique valley. (Note that regions giving tachyonic 
glueballs were assigned x^ — 00 for the purposes of Fig. ^). 

Physically one expects —00 as a ^ 0, since M should be forced to lie on the SU{N) 

group manifold M^M = 1 in the continuum limit. Since we have limited ourselves to the region 
> 0, we are some way from the continuum. It is remarkable, therefore, that such a strong 
signal for the Lorentz trajectory T is obtained. We can deduce the lattice spacing in units of cr on 
Ts if we assume no scaling violations; see Fig. ^. Naively one expects m to increase with a, but 
measurements of a vs. m are very sensitive to scaling violations. Fig. ^jis roughly consistent for the 
smallest m region, though scaling violations appear much stronger in 3 -|- 1 dimensions than they 
were in 2 + 1-dimensions ||l^ . 
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Figure 4: Charts of coupling space with Xmin ^ the height variable. 
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m 

Figure 5: The variation of the lattice spacing a with link-field mass m along 7^, assuming exact 
scaling holds (this only holds, even approximately, at small m). 
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Con Coff 


M 


0++ 0+ 


1.07 1.07 


3.50(24) 


2++ 2+ 
2- 
1± 
0+ 


0.86 0.87 
Im — 
0.91, Im 0.72, Im 
0.99 0.98 


4.97(68) 
6.07(?) 
5.62(?) 

4.97(18) 


1+- 1± 

0- 


0.88,1.06 0.94,1.00 
1.07 1.07 


5.37(36) 
5.74(45) 



Table 2: The glueball multiple! components at the lowest overall = Xo; showing masses from 
^ 2P+P"(P = 0) and c from the dispersion relation (p^). "Im" indicates < and "-" 
indicates the quantity has not been measured. The jTa ±1 states are exactly degenerate. 



4 Glueball data 

4.1 Masses 

In the limit C — s- oo, in addition to x'^-boost invariance, the transverse lattice has exact discrete 
symmetries of the group D4 [Q, together with charge- conjugation C. Reflections Vi and V2, about 
xi = and X2 — respectively, can be used to estimate parity V — 7^i7'2'P3 from free-field estimates 
of the light-front parity T'a : a;+ «-> . 90-degree rotations x^ x^ are exact and can be used to 
distinguish the angular momentum projections jTis = 0,±1,±2 from each other. 

In Fig. ^, we plot the scaling behaviour of the lightest glueball masses M. along 7^, labelled 
by IJTsI^^ and grouped into would-be Spin-Parity-Charge-Conjugation multiplets J^'^ . We have 
estimated error bars on each datum by combining well-understood (and small) DLCQ and Tamm- 
Dancoff extrapolation errors (~ 0.03A^) with an estimate of systematic finite-a errors coming from 
violations of Lorentz covariance, using the formula 

M{xl)^\c-lU.. (22) 

Xo 

|c— l|av. is the averaged deviation of the speed of light in transverse directions x, relative to direction 
x^. The datum of lowest overall = Xo occurs at m = 0.044. Formula (p^, which is supposed 
to represent roughly the uncertainty of the intercept Ai of the dispersion relation, works well for 



nearly covariant glueballs in 2 + 1-dimcnsional studies where the 'exact' mass is known |11, [l(j| . 
It is also consistent with the magnitude of scaling violations and Lorentz multiplet splittings in 
the present case. The various components of a Lorentz multiplet in Fig. ^ become rapidly more 
covariant, measured by their isotropy and degeneracy, as m is reduced]^ 

The final mass estimates given in the introduction are displayed in more detail in Table |^. They 
are taken from the point of lowest x^ = Xo = 13.36 on the scaling trajectory. We are confident that 

^ Naively, one might expect the colour-dielectric expansion to break down at the very smallest m, as the link fields 
become very light. However, the eneffry harrie r to production of more links comes principally from gauge self-energy 



contributions at these small masses ||?|, |32| . 
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0.1 0.2 0.3 0.4 





Figure 6: The variation of glueball masses with hnk-field mass m along Tg. Tiic open-symbol data 
for the 2++, which arc plotted only at = Xoi ^^111 too inaccurate for error estimates. Data at 
given m are displaced horizontally for clarity. 
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2 4 6 



Figure 7: Heavy-source potential Vqq vs. source separation R — [bx^Y + ((^a;^)^ + [bx^Y . Shown 
is a fit to the measured potential at zero transverse separations = Sx^ = 0. The 'off-axis' data 
are: (O) Sx^ = 1.22, Sx^ = 4.95, Sx'^ = 0; (X) Sx^ = 4.89, dx^ = 4.95, fa^ = 0; (+) Sx^ = 0, 
Sx^ = 4.95, Sx^ = 4.95. Units are {G^N)~^/^. 



the value of and error on these results is a reasonable estimate of the continuum result, provided 
the relevant observable exhibits a scaling region along 7^ and the error (|2^) is small. We emphasise 
that scaling alone is not sufficient for a good estimate of the continuum result. Table || gives some 
idea of where the remaining errors lie. In particular, the T'l = — 1 components of the tensor glueball 
are still dominated by lattice artifacts, requiring further terms in the Hamiltonian presumably. Our 
lightest ^ state lies at ~ 8^/a■, compared to much lower SU{3) Euclidean lattice Monte Carlo 



results 1 20 ; but we have no measure of its covariance and suspect a very large error on our result. 

It is interesting to plot glueball masses in pure gauge theory versus since this is supposed 

to be the relevant expansion parameter about N = oo |2^ ; see Fig. |l| in the introduction. There is 
remarkably little variation of glueball masses with the number of colours. It gives further support 
to the notion that = 3 is close to A'^ = oo in many situations. In particular, popular flux-tube 
and string models of the soft gluonic structure of hadrons are typically more appropriate to the 
large- iV limit of QCD (or are independent of N). Figure ^ indicates that these models should give 
worthwhile quantitative approximations to = 3 pure QCD . One can even go so far as to derive a 
formula for glueball masses at any N, by fitting those at = 2, 3, . . . , oo to polynomials in 
this might be a useful guide for higher- A' ELMC studies |^ . 

4.2 Covariance 

In searching for 71 we used data from the heavy-source potential and glueball dispersion near P = 
(the speed of light). At the — Xo point on 7^, the potential is plotted in Fig. |^. The potential 
between sources separated by distance K only in the x^ direction has been fit to 

= 0.084i?VG2]V-|-0.259- ^^^EL . (23) 
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2P+P- 



0.5- 




a|P| 



Figure 8: Dispersion relation of 0+"'" glueball. Straight line is exact relativistic form, dots are data 
at = Xo point on 7^ (black for direction P cx (1,0), gray for P cx (1,1)). The dashed line is a 
'lattice' dispersion cx ^1 ~ cosa|P| along P cx (1, 0), normalised to the correct slope at P = 0. 



The remaining points in Fig. are 'off-a;'^-axis' samples of the potential at a general separation 
i?, whose deviations from the value given by ( |2^ ) are used in the test. As can been seen, 
the rotational invariance of the potential is patchy, showing significant deviations at very small 
separations and/or along the (1,1) transverse direction (see discussion in Section pT^ ). 

Only the long range form of the glueball dispersion was used as a variable in the x^ test. As a 
further consistency check, we can also plot the dispersion throughout the Brillouin zone < |P| < 
7r/a. As the example of the lightest 0++ glueball illustrates (Fig. ^), improvement towards the full 
relativistic form seems to occur on 7^. 

4.3 Wavef unctions 

From the explicit glueball eigenfunctions one can extract various measurements of glueball struc- 
ture. While these results are still somewhat academic until one couples quarks to allow the glueballs 
to decay, they are still important. One can check whether the constituent structure of boundstates 
is valid. The probability distribution of transverse shapes is a useful measure of this (see Table ^). 
If eigenfunctions were to contain large mixtures of different numbers of link-partons, it would call 
into question the validity of the colour-dieletric expansion. 

Another interesting quantity is the distribution of longitudinal momentum P+ among the link 

^"The results in this subsection are from an earlier calculation that used a slightly different choice of test 
than that appropriate to Table hi The difference in the resulting wavefunctions is negligible for the purposes of our 
discussion in this subsection. 
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scalar 




tensor 




veetor 


Shape 


{0+} 


{2+ 


2- 


1± 


0+} 




0-} 


m 


0.019 


0.016 








0.014 


0.810 





h 


0.003 


0.518 








0.535 


0.039 





h 


0.002 


0.369 








0.361 


0.016 





h 


0.483 





0.882 


0.923 


0.017 


0.052 


0.049 


b 


0.400 





0.024 





0.003 


0.069 


0.860 


other 


0.093 


0.097 


0.094 


0.077 


0.070 


0.014 


0.091 



Table 3: Probability distribution of transverse shapes for various iJTal^^ components of the low-lying 
glueball eigenfunctions (following Table ||) . The shapes are denoted by the coupling constant for 
that shape in P^; for example, m denotes the TrjMrMj'} operator. Contributions < 0.001 are 
neglected. 

partons. In Fig. ^ we plot the quantity 

Gd(x,l/a) = I dx-e-'-P^-~{^{P+)\Tr{d-Mrd^M^^}\^{P+)) , (24) 

which measures (in A- = gauge) the probability of finding a link-parton carrying momentum 
fraction x of the glueball momentum P^. It depends upon the transverse normalisation scale 
through a, though our piece of the Lorentz trajectory is too short to reliably see physical evolution 
of Gd with scale, a 0.65 fm for the data shown. 

Gd is related to the gluon distribution; it becomes the gluon distribution in the limit a — > 0. 
Moreover, since M is some collective gluon excitation and the momentum sum rule is satisfied, one 
would naively expect the gluon distribution at a general scale a to be softer than Gd- In this case, 
the 0++ glueball (Fig. ^) does not resemble gluonium (two-gluon boundstate), which would have a 
distribution symmetric about x — 0.5. The fact that all components of the tensor's distribution look 
like the 0"*"^ is also interesting. It may indicate that the tensor's angular momentum comes mostly 
from gluon spin; however we caution that not all components of our tensor are yet behaving in a 
covariant fashion. Once light quarks are coupled to the problem, these distributions should have 
distinctive experimental signatures. 

5 Conclusions 

In this paper we have carried through the program, first suggested by Bardeen et al. [Q| , for solving 
the boundstate problem of pure QCD on a transverse lattice. Three of our developments were 
crucial to this success: inclusion of higher Fock sectors via DLCQ; inclusion of the full set of allowed 
operators occurring to lowest order of the colour-dielectric expansion of the light-front Hamiltonian; 
efficient semi-analytic methods for computing and extrapolating many-body light-front Hamiltonians 
p^ . In addition, we combined the analysis of the pure glue sector with an analysis of the heavy 
source sector, which is needed not only to provide an absolute scale but also to pin down the correct 
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Lorentz-covariant Hamiltonian in each sector. The upshot of all these developments is that we were 
able to clearly identify a scaling trajectory in the space of couplings on which Lorcntz covariance of 
observables is enhanced. 

We believe that this trajectory is an approximation to the exact Lorentz-covariant, scaling tra- 
jectory in the infinite dimensional space of all Hamiltonians. Strong evidence in favour of this comes 
from the scaling behaviour of glueball masses, whose values agree with continuum results (extrapo- 
lated to large- A^) from established methods. It is largely straightforward to relax the approximations 
we made so far, including the large- iV limit, to increase the scope and accuracy of our results. 

The light-front wavefunctions are entirely new results. Their accuracy is largely undetermined 
and, in their present form, they cannot easily be used to compute an experimentally measured 
quantity. However, their interesting structure might provide a source for phenomenological models. 
In fact, given the accuracy of the glueball masses and covariance of the dispersion relations, we 
would be very surprised to find that they were inaccurate. 

To make contact between light-front wavefunctions and the wealth of hadronic experimental 
data, the next step must introduce finite- mass propagating quarks. Although overviews of the 
standard Kogut-Susskind Q and Wilson fermions on a transverse lattice have been given, no 
first-principles calculations in the manner of this paper have been completed. By adding observables 
that test chiral symmetry (and parity) we can hope to find a similar scaling trajectory. Such 
calculations are now underway. 
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